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I. INTRODUCTION 



Dynamical chiral symmetry breaking and confinement are fundamental properties of 
QCD. In high energy processes such as deep inelastic scattering, quarks behave almost mas- 
selessly. However at low energies the observed hadron spectrum suggests that light quarks 
acquire large, dynamically generated masses through their interaction with the gauge sector 
of QCD. Quarks and gluons carry color charge and are not observed as asymptotic states, 
only occurring inside colorless bound states, the hadrons. The mechanism for such con- 
finement in QCD is still not understood and it is not known whether a gauge invariant 
formulation even exists. However, in the framework of a quantum theory, physical degrees 
of freedom are necessarily subject to a probabilistic interpretation implying unitarity and 
positivity; the physical part of the state space of QCD should be equipped with a positive 
(semi-) definite metric. Therefore one way to investigate whether a certain degree of free- 
dom is confined, is to search for positivity violations in the spectral representation of the 
corresponding propagator. Negative norm contributions to the spectral function signal the 
absence of asymptotic states from the physical part of the state space of QCD and are thus 
a sufficient (though not necessary) criterion for the confinement of the particle in question. 

Neither confinement nor dynamical chiral symmetry breaking can be accounted for at any 
finite order in perturbation theory. These phenomena can only be explored in genuinely non- 
perturbative approaches such as those provided by lattice Monte-Carlo simulations (see e.g. 
Ref. 13) and the Dyson-Schwinger, Green's functions approach (see e.g. Refs. jllUI^). 
Both approaches have their own strengths and weaknesses. Lattice simulations are the only 
ab initio calculations available so far. They contain the full non-perturbative structure of 
QCD but are limited by the enormous computational effort they require and by uncertain- 
ties in the infinite volume and continuum extrapolations that are needed to connect with 
the physical world. Furthermore, the implementation of small quark masses in most lattice 
simulations is computationally very expensive and, as yet, state-of-the-art calculations use 
light quark masses 6-10 times the physical values, thus necessitating a further extrapola- 
tion. On the other hand, the Dyson-Schwinger equations for the propagators of QCD are 
continuum-based and can be solved analytically in the infrared but must be truncated to 
obtain a closed, solvable system of equations jE H, 01- Recently, a concerted effort has 
been made to combine the strengths of these two approaches and quite definite statements 
on the infrared behavior of QCD have emerg ed [iiQIillii. In this work we will apply a 
similar strategy to explore the analytic structure of the propagators of QCD from solutions 
in the spacelike Euclidean momentum region. 

This paper is organized as follows: In Sec. |H] we briefiy review the connection between 
positivity and confinement and outline the method we will use to investigate the analytic 
structure of the propagator in the timelike momentum region. In the third section we 
investigate positivity violation in the gluon and quark propagators which are obtained as 
solutions of Dyson-Schwinger equations in the truncation scheme of Refs. jiol 11 1. We find 
clear evidence for positivity violations in the gluon propagator. The origin of these positivity 
violations is a branch point at = 0, followed by a cut along the real timelike axis. For 
the quark propagator we find no positivity violations as long as a certain non-perturbative 
Dirac structure is included in the quark-gluon vertex. This Dirac structure is dictated by the 
Ward-Takahashi identity in QED, and is also likely to exist in QCD because of the similar 
nature of the corresponding Slavnov-Taylor identity. In Sec. IIVI we seek parameterizations 



of the quark propagator. We investigate the ability of a number of meromorphic ansaetze 
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to reproduce lattice data for the quark propagator. All the fits share the property of either 
a dominant real pole or a pair of complex conjugate poles very close to the real momentum 
axis. We also show that one can reproduce both the Dyson-Schwinger solutions and the 
lattice data by various parameterizations with branch point singularities, rather than poles. 
We give a summary of our results in the last section. 

II. POSITIVITY AND CONFINEMENT 

One of the most intricate problems in quantum field theories is the separation of physical 
and unphysical degrees of freedom. In QCD this problem is directly connected with the 
issue of confinement, since we are searching for the mechanism which eliminates the colored 
degrees of freedom from the physical subspace, Vphys, of the state space of QCD. In order to 
ensure a probabilistic interpretation of the quantum theory, Vphys is required to be positive 
semi-definite, whereas the total state space of QCD in covariant gauges has an indefinite 
metric. 

A possible definition of a positive definite subspace, Vphys, is given in the framework of 
the Kugo-Ojima confinement scenario 0*]. Assuming the existence of a well-defined BRST 
charge operator, Qb, the space of physical states is defined by 

^phys = {\phys) : QB\phys) = 0}. (1) 

Given the assumption of a well-defined, i.e. unbroken, global color charge, Q"^, it has been 
shown that the physical state space Vphys only contains color singlets, i.e. {phys\Q'^\phys) = 
[l^ Esj ]. In Landau gauge this assumption, the Kugo-Ojima confinement criterion, can 
be translated into the requirement that the ghost propagator should diverge more strongly 
than a simple pole at zero momentum 

In this scenario, longitudinal gluons as well as ghosts are removed from the physical 
spectrum of QCD by the BRST quartet mechanism (see e.g. Ref. ^^). The colored states 
are BRST-quartet states, consisting of two parent and two daughter states of respectively 
opposite ghost numbers. The latter states are BRST-exact and thus BRST-closed (due to 
the nilpotency of the BRST transformation). The BRST daughters are orthogonal to all 
other states in the positive definite subspace and thus do not contribute to physical S'-matrix 
elements. The parent states belong to the indefinite metric part of the representation space 
and are thus expected to violate positivity. Members of the elementary quartet related to 
gauge fixing are the ghosts, the antighosts and longitudinal gluons. 

As the two parent states of a quartet belong to the indefinite metric part of the com- 
plete representation space, violation of positivity would provide evidence for the correctness 
of the Kugo-Ojima picture. E.g. positivity violation for transverse gluons indicates that 
transverse gluons are BRST-parent states with gluon-ghost states as daughters. The cor- 
responding parents of opposite ghost number are gluon-antighost states with a mixture of 
gluon-ghost-antighost and 2-gluon states as daughters. A similar construction for quarks 
would consider quarks as BRST-parent states with quark-ghost states as daughters, and cor- 
respondingly, quark-antighost states as second set of parents and a mixture of quark-ghost- 
antighost and quark-gluon states as second type of daughter states. Thus an investigation 
of (non-)positivity of transverse gluons and quarks allows us to understand in more detail 
confinement via the BRST quartet mechanism 

In order to complete the proof of confinement in this scenario one must still demonstrate 
the appearance of a mass gap in Vphys and the violation of cluster decomposition (see e.g. 
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Ref. 3; 121 and references therein) for colored states. Both requirements are related to the 
area law in the Wilson loop and, correspondingly, to a non-vanishing string tension in the 
quark- ant iquark potential. 

At this point we note that the basic assumption of the Kugo-Ojima confinement scenario 
still seems far from being proved: BRST-symmetry is a perturbative concept and it is not 
clear whether the symmetry remains unbroken in non-perturbative QCD jl8|. Furthermore, 
although clear evidence for a linearly rising potential between static quarks has been found in 
quenched lattice simulations (see Ref. jl^ and references therein), a mathematical proof of a 
violation of cluster decomposition is not at hand. Nonetheless, the Kugo-Ojima confinement 
criterion in its Landau gauge formulation has been tested in Dyson-Schwinger studies and 
in lattice simulations. Both methods agree very well even on a quantitative level and find a 
strongly diverging ghost propagator at small momenta [13, Ell HO, Hll, li^ . 

The Kugo-Ojima scenario is one particular mechanism that ensures the probabilistic 
interpretation of the quantum theory. However, even if it were eventually shown not to 
be appropriate, it is apparent that there is some mechanism which singles out a physical, 
positive semi-definite subspace in QCD. This suggests another criterion for confinement, 
namely violation of positivity. If a certain degree of freedom has negative norm contributions 
in its propagator, it cannot describe a physical asymptotic state, i.e. there is no Kallen- 
Lehmann spectral representation for its propagator. 

Within the framework of a Euclidean quantum field theory (which is used throughout 
this work) positivity is formulated in terms of the Osterwalder-Schrader axiom of reflection 
vositivitv Ual i . (F or a thorough mathematical formulation of the axiom the reader is referred 
to Refs. |24 liql). In the special case of a two-point correlation function, A{x — y), the 
condition of reflection positivity can be written as 

J d'^x d'y fix, -xo) A{x - y) f{y, yo) > , (2) 

where /(x, xq) is a complex valued test function with support in {{x, xq) : xq > 0}, i.e. for 
positive times. After a three-dimensional Fourier transformation, this condition implies 



oo 

I 



dtdt' f{t',p}A{-{t + t'),p)f{t,p} >0. (3) 



Provided there is a region around to = ~(^ + t') where A{tQ,p) < 0, one can easily find a 
real test function f{t) which peaks strongly at t and t' and thereby demonstrate positivity 
violation. For the special case p = 0, the Osterwalder-Schrader condition, Eq. Q, can be 
given in terms of the Schwinger function, A{t), defined by 



oo 

74, 



A{t) := j d'x I Al-/(^P^+^-P^a{p') = dp,cos{tp,)a{pl) > 0, (4) 



where cr(p^) is a scalar function extracted from the corresponding propagator. For the prop- 
agator of transverse gluons, cr(p^) is simply given by the renormalization function times the 
tree-level expression (see Eq. ()14|1 below) and we denote the corresponding Schwinger 
function by Ag{t). The quark propagator can be decomposed into a scalar and a vector part 

S{p) =: z^a,{p') + a,{p'), (5) 
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leaving us with two scalar functions, cyvip^) and crs(p^), to form two Schwinger functions, 
A„(t) and A^(t). 

Two simple examples for the analytic structure of a propagator in a quantum field theory 
are a real pole and a pair of complex conjugate poles. These highlight the paradigmatic 
behaviors of the Schwinger function, Eq. (0)). In the following, we always discuss the propa- 
gators and the functions (Js,v{p'^) in terms of the Lorentz invariant complex momentum, p^. 
Our notation is such that positive real values, > 0, correspond to spacelike momenta. 

(I) Real pole. The propagator of a real, massive, scalar particle has a single pole on the 
real timelike (p^ < 0) momentum axis. In this case the propagator function is given by 
cr(p^) = + m^) and it is easy to see from Eq. Q that the Schwinger function decays 
exponentially, 

A(t)~e-™*, (6) 

and is positive definite. For a bare propagator, the pole mass, m, is the same as the bare 
mass occurring in the Lagrangian. However, for an interacting particle, the pole mass can 
have both tree level and dynamically generated contributions. The real pole corresponds to 
the presence of a stable asymptotic state associated with this propagator. This does not 
imply that this state corresponds to an observable physical particle: provided the Kugo- 
Ojima scenario holds, all states belonging to a quartet representation of the BRST- algebra 
are excluded from the physical subspace, Vphys, which contains only colorless singlets. Thus 
two-point correlations of colored fields may develop real poles in momentum space without 
contradicting confinement [i^. In lattice calculations 0| and other non-perturbative ap- 
proaches j27|, the exponential decay in Eq. © is used to extract hadron masses and other 
observables from the large time behavior of appropriate correlators. 

(II) Complex conjugate poles. Another possible analytic structure for a propagator is 
a pair of complex conjugate poles with "masses" m = a ± ib. As has been discussed in 
detail in Refs. (28|, such a propagator could describe a short lived excitation which decays 
exponentially at large timelike distances. Furthermore, it has been argued [28.J that although 
causality is violated at the level of the propagators, the corresponding S-matrix remains 
both causal and unitary. Such complex conjugate poles lead to oscillatory behavior in the 
Schwinger function, A{t). Specifically, 

A{t) ^e-"' cos{bt + 6). (7) 

In this case one has negative norm contributions to the Schwinger function and the effective 
mass, 

^eff(t) = ^ (8) 

(defined in analogy to the real pole case, Eq. (0)) exhibits periodic singularities. Therefore 
the associated state (if there is any) must be an element of the unphysical subspace. Under 
the assumption of an unbroken BRST symmetry, this state must be a member of a BRST 
quartet, and the corresponding excitation is confined. 

Complex conjugate poles have been found for the fermion propagators of QED3 j29l |. 
QED4 (see e.g. H), and QCD |i3, H H H H Q in a variety of truncation schemes. 



In a number of these studies, the authors have discussed whether the observed positivity 
violations are genuine pro per ties of the theory related to confinement or artifacts of the 
truncation schemes (iollsi.lsGl 37|. As examined in the following section, it is our contention 



that dominant complex conjugate poles are indeed an artifact of the rainbow (bare vertex) 
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truncation of the quark Dyson-Schwinger equation and that, at least in Landau gauge, 
confinement through positivity violation in the quark propagator is not manifest. Complex 
conjugate propagators are also known to be practicable in light-cone dominated processes 
[ssj l and have recently been investigated in terms of the solution of the Bethe-Salpeter 
equation [s^. It has also been suggested that the gluon propagator may have such an 
analytic structure jH, 4^ 41, i^l- This possibility has been investigated in Refs. |43l. li^. 



Here, a note on positivity for the propagator of a Dirac field is in order. A dispersion 
relation representation of a fermion propagator in Minkowski space reads 







' — s + ze 



and positivity amounts to the requirements that for s > 

Pvis) > , and ^/s p^{s) - ps{s) > . (10) 

It is obvious that for a free Dirac field of mass m one has 

Pv{s) = 5{s — m^) , and Ps{s) = m S{s — m^) , (11) 

and thus ^/spy{s) — Ps{s) = 0. For an interacting Dirac field with physical asymptotic states 
and mass m one expects ps,v{s) = for s < m^. For s > m^, Eq. (jiup has to be satisfied. 
This requirement is automatically fulfilled if the stronger constraint 

mpr,{s) > ps{s) , (12) 

holds. 

Given the linearity of the different types of integral transforms relating <Js,v{p'^), Ps,v{s), 
and Ag^yit) to each other, one can conclude that cr^(p^) must be multiplied by a typical mass 
scale before being compared to as{p'^)- Thus, positivity violations can be signaled either in 
(j^(p^) alone, or in appropriate linear combinations of (Tvip^) and asip^)- We also consider 
the Schwinger function associated solely with o"s(p^), since it can be calculated with greater 
numerical accuracy. In general, oscillatory behavior in As(t) signals oscillatory behavior in 
At,(t) as well. 

Using the corresponding Schwinger functions, we can search for possible positivity vio- 
lations and investigate the analytic structure of the gluon and quark propagators of QCD. 
The t-dependencies of these Schwinger functions are determined by the analytic properties 
of the propagator, and, for large t, are dominated by the singularity closest to = 0. A 
complementary, direct method of determining the analytic structure is to solve the corre- 
sponding Dyson-Schwinger equation over a large region of the complex momentum plane. 
However, from a numerical point of view, such a procedure is very expensive and is not fea- 
sible with the resources currently available to us. Furthermore, there is good evidence from 
an investigation of QED3 that both methods agree very well jl^. We are thus confident 
that the Fourier transformation method is able to determine the qualitative behavior of the 
propagators. 

To complete this discussion we note that the conversion of a tree-level pole into an alge- 
braic branch point with exponent larger than one is also known for certain approximations 
to the fermion propagator of QED4 (see, e.g.^ supplement S4 in Ref. 0] and references 
therein). This type of singularity, {p^ + m^)^^^"/'^, is related to the soft photon cloud. The 
examples discussed in this section (real poles, complex conjugate poles, or branch cuts) 
will form the basis of our investigation of the analytic structure for the quark and gluon 
propagators. 
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III. SOLUTIONS OF THE PROPAGATOR DYSON SCHWINGER EQUATIONS 
OF LANDAU GAUGE QCD 



In this section we present solutions of the coupled set of Dyson-Schwinger equations 
(DSEs) for the ghost, gluon, and quark propagators in Landau gauge and investigate some 
of their analytic properties. In order to keep this paper self-contained, we first briefly review 
the DSE truncation scheme developed in Refs. [lO|, which is used to determine the 



propagators for Euclidean spacelike momenta, i.e. for real > 0. It is important to note 
that the behavior of the propagators for p"^ — > 0+ is extracted analytically. 

The DSEs for the quark, gluon and ghost propagators are derived from the QCD gen- 
erating functional with gluon field configurations restricted to the first Gribov region j40|. 
In a recent work it has been argued that such a prescription is sufficient to eliminate the 
effects of Gribov copies from correlation functions |4^. Furthermore, the DSEs are not 
affected by imposing such a boundary condition on the generating functional of the gauge 
fixed theory because the Gribov horizon is a nodal surface for the integrand of this func- 
tional integral. Instead, the ghost two-point function has to satisfy the so-called horizon 
condition 9], i.e. the ghost propagator has to diverge more strongly than a simple pole for 

^ 0^. This condition (which in Landau gauge is formally equivalent to the Kugo-Ojima 
confinement criterion discussed in the preceding section) turns out to be enforced by the 



ghost DSE [8|, |47|, |48[ and is thus fulfilled by the DSE solutions in the truncation scheme 
that we employ. 

A graphical representation of the DSEs for the ghost, gluon, and quark propagators is 
given in Fig. Q and their full form can be found in Ref . jsf . In Landau gauge (which is used 
throughout this work), the renormalized ghost, gluon and quark propagators, Dg(p, /i), 
Dfj,v{p,lj), and S{p,fi), respectively, are given in terms of scalar functions by 

DgVP^N = 5 ' 13) 

^/..(p, = (^V --^) — j2 — ' (14) 
^(P'/^) = =■■ ^M^^^') + -s{p^^^'). (is) 

All these propagators are diagonal in their respective representations of SU(Ai'c), so their 
color structure has been suppressed for simplicity. The dependence on the renormalization 
scale, fi, is given explicitly for later use. Here, G{p'^, /i^) and Z^p"^, ij?) are the ghost and gluon 
dressing functions, respectively, and A{p'^,fi'^) and B{p'^,ix^) are the vector and the scalar 
parts of the inverse of the quark propagator. The functions most relevant for our study 
of positivity are Z(p^, /i^)/p^, (Ts(p^,/i^) and (j„(p^,/i^). Note that the ghost propagator 
trivially violates refiection positivity because of the way ghosts are introduced in Faddeev- 
Popov quantization [49]. 

Two renormalization-scale-independent combinations built from the scalar functions rep- 
resenting the different propagators are important for further discussion. First, M(p^) = 
i?(j9^, yU^)/A(p^, yU^) denotes the renormalization-point-independent quark mass function. 
Second, as has been demonstrated in Ref. 0|, a non-perturbative definition of the running 
coupling, is possible due to the non-renormalization of the ghost-gluon vertex in Landau 
gauge [50|. This results in the relation 

a{p') = a{^^^)G\p\^J?)Z{p\^^'). (16) 
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FIG. 1: Diagrammatic representation of the Dyson-Schwinger equations for the gluon, ghost, 
and quark propagators. The wiggly, dashed, and sohd Hnes represent the propagation of gluons, 
ghosts, and quarks, respectively. A filled blob represents a full propagator and a circle indicates a 
one-particle irreducible vertex. 

In the following we investigate the full (unquenched) system of DSEs and also the quenched 
approximation to them in which quark loops are neglected, removing the back-reaction of 
the quarks on the ghost and gluon system. 



8 



A. Truncation scheme 



Both the quenched and the unquenched system of ghost, gluon, and quark DSEs have 
been solved numerically in Refs. 0, in a truncation scheme which neglects the effects 
of the four-gluon interaction and employs ansaetze for the ghost-gluon and the three-gluon 
vertices such that two important constraints are fulfilled: the running coupling, a{p^), is 
independent of the renormalization point and the anomalous dimensions of the ghost and 
gluon propagators are reproduced at the one-loop level for large momenta. In order to study 
the effects of violating gauge invariance by these truncation assumptions, the gluon DSE 
has been contracted with the one-parameter family of tensors 

Vii\p) = 6,.-C^-f'. (17) 

In Landau gauge, a violation of gauge invariance manifests itself in the appearance of spu- 
rious longitudinal terms in the gluon equation, which in turn introduces dependence of the 
ghost and gluon dressing functions on the parameter (. The influence of these longitudi- 
nal terms has been examined in Ref. by varying ^ and found to be surprisingly small. 
Further technical details of the truncation scheme in the Yang-Mills sector are relegated to 
Appendix A where we also discuss the dependence of our analysis on these details (see also 
Refs. 

Employing asymptotic expansions for the propagators at small momenta, the untruncated 
ghost and gluon DSEs can be solved analytically for 0+ ji^]. One finds simple power 
laws, with exponents related as 

Zip^f^') ~ {pVf^'r, (18) 
G(/,/i^) ~ {pVf^Y\ (19) 

for the gluon and ghost dressing functions. The value of the exponent n depends somewhat 
on the details of the employed truncation scheme. In certain truncations it can be calculated 
analytically and it will depend on the parameter ( [l^. The tensor vjxi~^^ projects onto 
the purely transverse part of the gluon equation, and in this case the solution k, = (93 — 
v/T20T)/98 ^ 0.595 has been found in Refs. By varying 1 < ( < 4, infrared solutions 

with exponents in the range 0.5 < K < (93 - Vl201)/98 have been shown to connect to 
numerical solutions for all momenta |lOl |. A recent infrared analysis of the ghost and gluon 
DSEs employing the most general ansatz for the ghost-gluon vertex suggests the exponent 
n is in the range 0.5 < k < 1 (which is further restricted to 0.5 < k < 0.7 after 
constraints on the value of the running coupling are taken into account). A first attempt 
to include the two-loop diagrams in the gluon DSE also results in very similar values for 
the infrared exponent jHH and in Ref. |4y| it has been shown that the two-loop diagrams 
have no effect on k. Finally, exact renormalization group equations have recently been 
employed in a complementary investigation [s^ of the infrared behavior of the gluon and 
ghost propagators with a resulting value for k in agreement with those above. These varied 
investigations all indicate that the Landau gauge gluon propagator vanishes as p"^ — >• 0"*" and 
predict an exponent 0.5 < k < 0.7. 

For the subsequent discussion, it is important to note that the exponent k is very likely 
an irrational number. The relation of the exponents in Eqs. (fTHj) and (fT^ results in an 
infrared finite strong coupling independent of the value of k, c.f. Eq. (fT^ . For transverse 
projection, the value is given by a(0) = 8.915/iVc. 



9 



The DSE for the quark propagator S{p, n) is given by 



(20) 

where Z2 and Zip are the quark wave function- and quark-gluon vertex-renormahzation 
constants, respectively, and J represents a translationally-invariant regularization charac- 
terized by a scale, A. The momentum routing is = q—p, and the factor Cp = {N^ — l)/2Nc 
stems from the color trace of the loop. 

In addition to the quark and gluon propagators, Eq. (pUj) involves the quark-gluon vertex, 
Ty{q^p] jj). This vertex is, in principle, determined by its own DSE jSSf involving various 
(n < 5)-point correlators. However, the solution of such higher-order DSEs is difficult even 
in the simplest situations jH^I and we avoid the problem by making an ansatz for T^{q,p; /i). 
As the structure of this vertex turns out to be crucial in our analysis of positivity violations 
in the quark propagator, we explore its construction in some detail. 

A reasonable ansatz for the quark-gluon vertex has to satisfy at least two constraints: it 
should guarantee the multiplicative renormalizability of the quark propagator in the quark 
DSE, and it should at least approximately satisfy its non-Abelian Slavnov-Taylor identity. 
It has been shown in Ref. [ill that the construction 



r,(g,p;/i) = V:'^'{q,p-fi)W-^'^%p\q',k'-f,) , (21) 

with 

W--'^\p',q^e;fi) = G'ie,f^') Z3(/x^A2) , (22) 

V:''Kq,p;fi) = T^''{q,p;f^) 

A(p2,/i2) + A(g2,/i2) ^ B{p^fi^)-B{q^fi'^ 



Iv + ^ 2 5 + 

{i> +i){p + q)u 



2 p^-q^ 
A(p^/.2)-^(g^/i2^ 



2(p2 - g2) 

X P' + g' (23) 

{p^ -q^y + {M^{p^) + M^{q^)y' ^ ' 

and Z3 being the ghost wave function renormalization constant, satisfies these requirements. 
Here it is assumed that the non-Abelian part of the vertex, W^"'^'^''{p'^,q'^,k'^; n), can be 
factored out from the Dirac structure, and that the Dirac structure is given by F^^fq, jo; u\ 
the Curtis-Pennington (CP) construction of the fermion-photon vertex in QED4 55|, |56[. 
Note that the dressing of the longitudinal part of the CP vertex is dictated by the Abelian 
Ward identity 

-zk, T'^^''''{q,p;fi) = S-\p,fi) - S-\q,fi) , (24) 

which results, among other things, in the appearance of a quark-gluon coupling term pro- 
portional to the sum of the incoming and outgoing quark momenta, 

p2. _ qz 
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Such a coupling, being effectively scalar, may at first sight appear to violate chiral symme- 
try, as, in contrast to the perturbatively dominant vector coupling proportional to 7,^, the 
expression ()25p commutes with 75. However, this scalar term only appears if chiral symme- 
try is already dynamically broken and is thus consistent with the chiral Ward identities. Its 
existence provides significant additional (self-consistent) enhancement of dynamical chiral 
symmetry breaking. Such a scalar coupling also appears in vertices that occur in system- 



atic improvements on the rainbow (bare vertex) truncation [57|, |58|, |59| . This term will be 
important in our investigations of positivity below. 

For comparison, we also employ a construction with a bare Abelian part of the vertex 
given by 

Vf^\p,q-li) = Z2(/i,A)7.. (26) 

In both cases the input from the Yang-Mills sector of the theory, i.e. the factors from the 
dressed gluon propagator and the non- Abelian vertex dressing can be combined to give the 
running coupling a{k^) = g^{fi) G^{k^,fj,'^) Z {k^ , fi"^) / Att according to Eq. (fTBj). Thus we 
arrive at the truncated quark DSE 

S~\p,f.) = Z,{^^')S,\p) + ^^ I d'q^(^6,.-'^y,S{q,f,)Vnq,p-f,). 

(27) 

In the quenched and unquenched calculations of the quark propagator we take a{k^) directly 
from the ghost and gluon equations. 



We also consider the solutions of the quark DSE in the model calculations of Refs. |60ll61 
6^ . There, only the leading 7^-part of the quark gluon vertex has been employed and the 
combination of the gluon and vertex dressing needed in the quark DSE has been modeled 
phenomenologically. With 7^ = 12/(llA''c — 2iV^ being the anomalous dimension of the 



quark propagator, we follow the authors of Ref. [61| and use the model 

^(q'^) =^nn^ F--?'/-' I 7r7^ [1 -exp(-gVm^2)] 
q^ u;^ ^ +g4ln[e2-l + (l + gVA^^jr 

with Aqcd = 0.234 GeV in the MS-scheme, Nf = 4 and the parameters rrit = 1.0 GeV, 
uj = 0.3 GeV, and D = 0.781 GeV^ fixed by fitting the chiral condensate and pion decay 
constant. Omitting the perturbative logarithmic tail, we also compare with the model of 
Ref. using a purely Gaussian interaction 

^ = >g^e-^/-% (29) 
q^ uj'^ 

with u = 0.5 GeV and D = 1 GeV^ 

Despite the fact that these models for the effective interaction were designed to be used 
in combination with a bare vertex, we also use them in conjunction with the CP vertex, 
T'^^. By comparing quark propagators that result from employing either direct input from 
the ghost and gluon sector or the model forms, Eqs. (j28p and (|29p . we are in a position 
to test whether the analytic properties of the quark propagator are more sensitive to the 
global strength of the quark-gluon interaction, to the overall shape of the (effective) running 
coupling, or to the details of the tensor structure of the quark-gluon vertex. First however, 
we will discuss the results of the numerical calculations for the gluon propagator. 
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B. Results for the gluon propagator for Euclidean momenta 
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FIG. 2: The solutionis of the quenched {Nj = 0) and unquenched {Nf = 3) coupled DSEs for 
the gluon dressing function, Z^p"^), are shown. The unquenched case with three massless flavors is 
compared to different fits (see text for details of the fits) . In the left panel these are displayed on 
logarithmic scales, in the right panel, on linear scales. Results from quenched lattice calculations 
[631 are given in the right panel. 

In Fig. 121 we display the numerical results for the gluon dressing function calculated with 
zero (quenched) or three (unquenched) flavors of massless quarks and transverse projection, 
C = 1 (c./. Eq. (fT7|) ). taken from Ref. ^3-^ diagram on the right of Fig. El the 

DSE results are compared to results from quenched lattice Monte-Carlo simulations jo^. 
The quenched DSE results are seen to agree well with the lattice data. In contrast, the 
unquenched DSE gluon propagator is significantly suppressed in the intermediate momentum 
region where the screening effects of quark-antiquark pairs become important. For both 
Nf = and 3, there are two qualitative properties that we can extract from these results: 
the analytically calculated infrared behavior given by Eq. (fTH|) . and a maximum around 
~ 1 GeV, followed by relatively flat momentum dependence above this scale. 

The behavior of the gluon dressing function in the infrared is captured by either of the 
irrational functions^ 



which are exact in the infrared limit (c./. Eq. ()18p) and which play a role when it comes 
to the interpretation of our results for the gluon Schwinger function, Ag(t). The value 



^ As can be inferred from Refs. |lCtllll |. changing the projection of the gluon equation in the range 1 < C < 4 

leads only to quantitative changes in the gluon and ghost renormalization functions. 
^ From here on we shall suppress the renormahzation scale dependence (whenever possible) for concision. 




(30) 



(31) 
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for the exponent k ~ 0.595 in these fits is taken from the infrared analysis of the DSEs. 
Note that for k — 1, the form of Eq. ()30|) becomes identical to the Gribov form proposed in 
Refs. ^1,16^]. The normalization parameters w/, wn and scales A/, A// are chosen such that 
the Schwinger function of the (C = 1, A^/ = 3) numerical gluon propagator is reproduced 
by the Fourier transforms of the fits (the value of these parameters are given below). Our 
fits with these irrational functions Z^jjip'^) are shown in Fig. |21and clearly reproduce the 
behavior of the DSE gluon propagator for very small momenta but deviate significantly from 
the dressing functions at momenta above ~ 400 MeV. 

To describe the behavior for larger momenta, we multiply the functions ZfnijP') by a 
function incorporating the known ultraviolet behavior. To this end we note that in Ref. 11 
the numerical running coupling has been fitted by^ 

2 ^ a(0) 47r f 1 Aqcd ^ . . 

"''^^ ^ 1 + pV^Icd /5o + P^lcD \HpV^Icd) p' - Kcd) ' ^ ^ 

In this expression the Landau pole has been subtracted as has been suggested in the frame- 
work of analytic perturbation theory ;65]. The value a(0) = 8.915/A^c is known from the 
infrared analysis and = (llA'c — 2Nf)/3. Using a MOM scheme and fitting only the 
ultraviolet behavior, a value Aqcd = 0.71 GeV has been given in Ref. jlll |. 
Identifying Ajji = Aqcd for simplicity, we utilize the fits 

ZiMp') = Zlnip") «fi7(p') ' (33) 
for further investigations, using the one-loop value of the gluon anomalous dimension, 7 = 
(— 13A'c-|-4A^/)/(22A'c — 4^"/). The quality of these fits can be seen in Fig.|21 For a discussion 
of the parameters used, see below. 

Employing a numerical Fourier transform routine, we can now calculate the Schwinger 
function, Ag{t) (defined by Eq. for the numerical solutions of the gluon DSE and for 
the various fits. The absolute values of the numerical Schwinger functions for Nf = 0, 3 
(using transverse projection C = 1) are displayed in Fig. |21 The spikes mark the time scales 
at which the Schwinger functions cross zero and negative norm contributions appear in each 
gluon propagator. One notes that the Schwinger function in the quenched approximation 
differs visibly from that for three flavors, despite the similarity of the corresponding gluon 
dressing functions for Euclidean momenta |llj. In particular, the typical time scale, marked 
by the zero of the Schwinger function, decreases from 5.2 GeV~^ to 4.4 GeV~^. We have 
also explicitly checked that different choices for the projection of the gluon equation and 
other minor details of the truncation scheme lead only to minor quantitative alterations 
(see Appendix A). All gluon Schwinger functions we have calculated from the results of 
the coupled DSEs show the same qualitative behavior, thus demonstrating that neither the 
details of the projection in the gluon equation nor the feedback of (a small number of) 
dynamical quarks^ have any significant influence on the overall analytic structure of the 
gluon propagator. We clearly observe positivity violations in the gluon propagator. This is 
the first major result of this work. 



^ In Ref. 11] two additional parameters a and b were used with a — 1.020 and b — 1.052. As the deviations 
from unity are completely insignificant we have fixed a = 6 = 1 here. 

The infrared (p^ O"*") behavior of the Yang-Mills sector of QCD is unaffected by the appearance of 
chiral quarks as long as the number of flavors is small enough to be in the confining and chiral symmetry 
breaking phase of QCD [ll|. 
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FIG. 3: The results for the absolute value of the gluon Schwinger function, Ag{t), corresponding 
to our numerical results from the DSEs are shown and compared to the fits in the infrared (left 
panel) and the overall fits (right panel). The spikes mark the time scales where the Schwinger 
functions cross zero and negative norm contributions appear. 



C. Analytic structure of the gluon propagator 



In the following we aim at an interpretation of our results in terms of the analytic structure 
of the gluon propagator in the timelike momentum region. As a first step we demonstrate 
that the infrared behavior of the gluon propagator, i.e. the behavior for — > O"*", is 
responsible for the non-trivial analytic structure. To this end, in the left-hand side of Fig. El 
the numerical results for the gluon Schwinger function (with C = I5 = 3) are compared 
to the infrared fits, Eqs. and (jHH). The fitted parameters are wj = 2.5, A/ = 400 
MeV for IR-fit I, Eq. and wu = 2.7, An = 420 MeV for IR-fit II, Eq. dSH). As we 

observed earlier, the fits only agree with the numerical gluon dressing function in the infrared 
momentum region. Nevertheless, in Fig. 01 we see that the agreement of the numerical 
Schwinger function with the Fourier transforms of each of these fits is excellent. It appears 
that the details of the intermediate and large momentum behavior of the gluon propagator 
have little or no influence on the qualitative analytical structure of the propagator in the 
"near-by" timelike momentum regime. In particular, the change in curvature at the bump of 
the gluon dressing function at a scale of ~ 1 GeV is not an important feature in this regard. 
In fact the crucial property of the gluon propagator is that it goes to zero for vanishing 
momentum. This can be seen easily as the relation, 

0.^(,.0)^ (34) 

(with D{p) = Z{jP')Ip^) implies that the propagator function in coordinate space, D{x), 
must contain positive as well as negative norm contributions, with equal integrated strengths. 

For fit I (Eq. (jHSI)) we have used two parameter sets, wi = 2.4, Aqcd = 500 MeV and 
wi = 2.0, Aqcd = 470 MeV. The first parameter set fits the gluon renormalization function 
better (especially in the ultraviolet) and the second set is optimized to fit the Schwinger 
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function. For fit II (Eq. ((221)) with the parameters wu = 2.5 and Aqcd = 510 MeV both 
the gluon renormahzation function and the Schwinger function are fitted very well. As the 
infrared fits I and II already reproduce the gluon Schwinger function it is no surprise that 
the complete fits, Eq. do even better, see the right-hand side of Fig. El As already 

stated, for the sake of simplicity we have used only one common scale, Aqcd, for the infrared 
and ultraviolet behavior. 

We are now in a position to deduce the possible analytic structure of the gluon propagator. 
We first observe that because of the infrared singularity of the ghost propagator, we expect 
a cut on the timelike momentum axis coming from the ghost-loop contribution to transverse 
gluons. As the ghost loop is the infrared dominant contribution in the gluon equation and 
therefore determines the infrared behavior of the gluon propagator, it is instructive to discuss 
the infrared fits to the gluon propagator first. The infrared fit I (Eq. (jHOI)) contains a branch 
cut on the negative axis while the denominator contributes a pair of complex conjugate 
singularities at 

p2 = ^2 g±W(2K) _ (35) 

The discontinuity across the negative axis is easily calculated. Writing p\ = (— p ± ie)A'j 
one obtains 

hm {D.ipl) - D:{p- )] = ^ M'^nn) p-''-' ^^^^ 

e-o^ JVF-^; 1 + 2 p2« cos(27r/t) + p4« ' ^ > 

with DiijP') = ZY{p'^)/p'^. This discontinuity rises from zero at p = to a maximum at the 
area of the pole locations and then rapidly decays as p becomes larger. 

In the infrared fit II (Eq. (jS^) the numerator and the denominator conspire to produce 
one cut^ over p^ G (0, — A^^). For the discontinuity we have (now for p\ = (— p ± ie)Ajj) 

r /n r n r 2U -2z^// sin(27rK) p^'^"^ 

\im{Dnip+) - Dnip^)} = -j^^ _ ^ (37) 

with Diiij?) = Zfj{p'^)/p'^ and for < p < 1 only. This rapidly diverges as | —A'jj (i.e. 
p — > 1) and then drops discontinuously to zero: there is no discontinuity for < — A^^. 

Whereas the location of the singularity p'^ = —Ajj in the infrared fit II is independent of 
the value of the exponent k, the location of the complex conjugate singularities of IR-fit I as 
well as the magnitudes of the cuts in both fits depend on k and therefore on the truncation 
scheme. Although the exact value of k depends on the details of the truncation, various 
methods suggest that the exponent n is in the range 0.5 < k < 1 [1, 0,E3, exactly 
this range which corresponds to the pair of complex conjugate singularities in IR-fit I being 
located on the first Riemann sheet in the left half of the complex p^-plane. In the limiting 
case K = 0.5, one obtains one real pole on the negative p^-axis in both fits, and in the other 
limit, K = 1, IR-fit I corresponds to a pair of purely imaginary poles, i.e. exactly the form 
proposed in Refs. 



To discuss the analytic structure of the full fits, Eq. (jHSj), we must also look at the analytic 
properties of the expression for the running coupling, Eq. ()32|). The Landau pole at spacelike 



^ Note that we have decided to take the ratio first and then we raise it to a non-integer power. Having 
this non-integer for the numerator and the denominator separately would lead to two overlapping branch 
cuts. However, we consider this an unnecessary complication. 
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= ^^QCD has been subtracted, so expression only has singularities on the timelike real 
axis. The logarithm produces a cut on this half-axis, and the corresponding discontinuity 
vanishes for 0~, diverges at = — Aq^^ and goes to zero for oo. In the fits 

I and II, Eq. (jSSl), the running coupling (Eq. ((221)) is raised to a non-integer power and 
multiplied by the infrared fits (Eqs. (j3U|) and (j31|) ). Thus, fit I also has a pair of complex 
conjugate singularities, at the same locations as those in Eq. (jHUI) . On the other hand, fit 

II has no non-analyticities other than the cut on the negative real axis. The discontinuity 
corresponding to the cut of the combination of the different factors in fit II is always positive, 
vanishes for 0", diverges at = — Ag^^ to +00 and falls to zero for —00. 

It is interesting to note the scale at which positivity violations occur. From Fig. El we 
determine that the zero crossing appears at t ^ 5 GeV~^ ~ 1 fm. This is roughly the 
size of a hadron and therefore the correct scale at which gluon screening should occur. 
One might speculate whether this represents an inherent, gauge invariant scale (as the 
locations of propagator poles are protected by Nielsen identities |66[), which is generated in 
the renormalization process. The pure power law Z{p'^) = (p^)^'^, which solves the system 
of DSEs in the case where the renormalization point /i is shifted to asymptotic values, is 
in perfect agreement with the scale-invariance of the underlying theory, corresponding to 
an infinite mass gap. Thus it is obvious that we can deduce the existence of a cut from 
the pure power laws, but we can not extract the related scale. This scale emerges from an 
interplay of infrared and ultraviolet properties of the theory, i. e. the transition of the gluon 
propagator from the infrared power law to its perturbative ultraviolet behavior. 

Before concluding this subsection we comment on what lattice Monte-Carlo simulations 
say about positivity violation in the gauge boson propagator. For unquenched QCD, nothing 
is known because the gluon propagator has not yet been calculated with dynamical fermions. 
The pure Yang-Mills gauge propagator has been calculated on the lattice for almost twe nty 
years following the pioneering work of Mandula and Ogilvie jo^], see e.g. Refs. im, Hill Eg 
and references therein. However, explicit observations of positivity violation have been 
elusive as statistical errors and finite volume artefacts cloud the issue. Nevertheless, many 
hints of negative norm contributions in the gluon propagator have been reviewed in j69|. 
Clear measurements of positivity violation have been made for the case of SU{2) and 
for the gluon propagator in three-dimensional Yang-Mills theory [t^ ■ 

Summarizing: the Landau gauge gluon propagator, as it results from the solution of cou- 
pled DSEs, displays positivity violations. This is in accordance with gluons being confined. 
The infrared behavior of the gluon propagator is analytically determined to be a power 
law. It has been demonstrated in Ref. p| that this behavior is stable under a broad range 
of possible dressirigs of the ghost-gluon vertex. Furthermore, strong arguments have been 
presented in Ref. |46[ for the existence of power laws in generalized truncations that include 
the four-gluon interaction. The power law behavior at small Euclidean momenta induces a 
cut on the real negative p^-axis, as can be seen clearly from our infrared fits. It is this cut 
which causes the observed pattern of positivity violation. Fitting the gluon propagator for 
all Euclidean momenta and the corresponding Schwinger function we are able to describe 
the gluon propagator with fit II, Eq. (j33|) . which has no singularities in the complex j9^-plane 
except for a cut on the negative real axis. 

Note that this fit contains essentially two parameters: the overall magnitude which, 
because of renormalization properties, is arbitrary,^ and the scale Aqcd- The infrared 



^ i.e. it is determined via the choice of the renormahzation scale fi and the normaUzation condition 
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FIG. 4: The quark mass function, M(p^), and the wave function renormahzation, Z^^p^), from 
quenched {Nf = 0) and unquenched {Nf = 3 chiral quarks) DSEs JA^. Results for the generalized 
CP vertex, Eq. (|23|) . and the bare vertex construction, Eq. (|26|) . are compared "with quenched 
lattice data in the overlap [z^l and Asqtad [l^l formulati 



tions 



exponent, k, and the anomalous dimension of the gluon, 7, are not free parameters: k 
is determined from the infrared properties of the DSEs and the one-loop value is used for 7. 
Therefore, we have found a parameterization of the gluon propagator which has effectively 
only one physical parameter, the scale Aqcd- Combined with the relatively simple analytic 
structure of fit II, Eq. (j^Sl), this gives us confidence that we have succeeded in uncovering 
the most important features of the Landau gauge gluon propagator. 



D. Results for the quark propagator 

In Fig. m we display the mass function, M(p^) = B{p'^)/A{p'^), and the wave function 
renormahzation, {p^) = l/y4(p^) (note the superscript / which differentiates this function 
from the gluon dressing function), of the quark pro pagator in the chiral limit, obtained 
from the coupled quark, ghost, and gluon DSEs We show quenched {Nf = 0) and 

unquenched {Nf = 3) results employing the generalized CP vertex, Eqs. (j2H)-(j2ni)- We 
also display the same functions calculated in the quenched approximation with the bare 
Abelian part of the quark gluon vertex, Eq. (j^Uj) . On the Euclidean real axis, both vertex 
constructions lead to qualitatively similar but quantitatively quite different results. The 
bare vertex approximation does not give enough chiral symmetry breaking and is clearly 
disfavored by recent quenched lattice data |72l. 173^ (also shown in Fig. H)). On the other 
hand, the results for the more elaborate vertex construction are well within the region 
suggested by the lattice calculations. 

The quantitative difference between the DSE solutions using the bare vertex and the CP 
vertex turns into a qualitative difference for the corresponding Schwinger functions. The 
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FIG. 5: The left diagram displays the absolute value of As(t) employing the bare vertex construc- 
tion in the quark DSE. The spikes correspond to zero crossings of the Schwinger function. These 
are absent in the diagram on the right where the results with the full CP vertex, Eq. (|2H) . are 
considered. The chiral limit results are shown for Nf = 3 and Nf = 0, together with the fits to the 
Schwinger function of the quenched DSE solution. Furthermore we compare to a calculation with 
only the two most important terms of the quark-gluon vertex. 



Fourier transformed scalar parts of the different quark propagators, As{t), are shown in 
Fig. El Similar results are obtained for the vector parts of the propagators, A„(t), though 
they are numerically less accurate.^ As in the case of the gluon propagator, we plot the 
absolute values of the Schwinger functions on a logarithmic scale. The results in the left 
diagram are obtained employing the bare Abelian part of the vertex, Eq. ()26|). Clearly these 
solutions exhibit the oscillatory behavior of Eq. ((7j), which is characteristic for a propagator 
with a pair of complex conjugate "mass-like" singularities. A fit of the expression in Eq. ^ 
to our result gives the locations of these singularities as msing = (209 + lOli) MeV. 

A completely different picture is obtained from the Schwinger functions constructed using 
the CP vertex, Eq. (j23|) . as can be seen in the diagram on the right of Fig. El Again we 
display results for the quenched case, Nf = 0, and the case of Nf = 3 chiral quarks. For 
Nf = we also make use of a fit to the running coupling as described in detail in Ref . ; 
for all practical purposes the results are almost indistinguishable from those obtained with 
the numerical a(g^) as a solution of the ghost-gluon DSEs. We find no traces of negative 
norm contributions, and in all fit of the oscillatory form of Eq. ((7|) to our results 

indicates that there is a singularity (almost) on the real timelike axis, with an imaginary 
part of at most 8% of its real part. The best fit is obtained for a real mass singularity at 
'^sing = 0.50 GeV. For both the bare and the CP vertex, the deviation of the fits from the 
data at small time scales suggests that there is additional structure in the DSE solution 
which the simple pole fits (Eqs. © and (|7|)) do not capture. We shall investigate this in 



In the chiral limit, the scalar part of the propagator, crs{p^), falls off like up to logarithmic corrections, 
because the function B{p^) falls off like f/p^, whereas <Ty{p^) falls off like This makes the Fourier 

transform of the scalar part easier to calculate numerically. 
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bare vertex 


ZjA^ + Ai3^j-term 


CF vertex 


YM a(A;^), unquenched, Nf = 3 


0.21(1) ± 0.10(1) i 


0.48(3) 


0.50(3) 


YM a(fc^), quenched [Nf = 0) 


0.21(1) ± 0.10(1) I 


0.48(3) 


0.50(3) 


fit A of Ref. [11], quenched 


0.209(4) ± 0.101(2) i 


0.48(3) 


0.50(3) 


fit B of Ref. [11] , quenched 


0.160(4) ± 0.076(2) i 


0.42(3) 


0.42(3) 


Maris-Tandy model [61], Eq. ^ 


0.55(1) ± 0.321(6) i 


0.96(6) 


1.1(1) 


Gaussian model [62], Eq. (|^ 


0.53(1) ± 0.167(3) i 


0.83(4) 


0.83(6) 


quenched QED (in units of lO^^A) 


1.79(6) ± 0.43(2) i 




1.51(9) 



TABLE I: Results for the fermion pole masses in the chiral limit for different interactions, as 
extracted from the behavior of the corresponding Schwinger functions. The quark masses are 
given in GeV, the QED4 results are given in units of the UV cutoff A, and are obtained with 
a = 1.2 for the bare vertex and a = 1.06 for the CP vertex. The errors are estimates of the total 
numerical error; the numerical error in case of a real mass singularity is dominated by the fact 
that, on a logarithmic scale, the Schwinger functions are not perfect straight lines. 



SecHVl 

By turning the different contributions in the vertex construction of Eq. (j^Hj) on and off, 
we have identified the term which is responsible for the qualitative differences between the 
left and right diagram of Fig. El In addition to the (dominant) vector part of the vertex 

:= ^ 7/. , (38) 

the presence of the scalar coupling AB^, Eq. (f^K|) . in the quark-gluon vertex is crucial for 
the substantial change in the analytic structure of the quark propagator compared to the 
truncation keeping only the vector part. Such a scalar term introduces additional feedback 
in the scalar self-energy, and its presence considerably enhances the amount of dynamical 
chiral symmetry breaking generated in the quark DSE. By varying the strength of this term 
compared to the leading SA^-piece of the vertex, we find that a reduction of this term 
by about 20% is enough to generate again positivity violations corresponding to dominant 
complex conjugate singularities. 

The question of positivity violation does not depend on the details of the input from the 
Yang-Mills sector of QCD. We obtain quantitatively similar results for the unquenched case 
with Nf = 3 chiral quarks, for the quenched approximation with the running coupling taken 
directly from the Yang-Mills DSEs and for different models for the running coupling 
As a check, we also employ the model interactions given in Eqs. ()28|) and (j^^ . Again we 
obtain evidence for a pair of complex conjugate singularities when a bare vertex is used and 
a singularity on the real timelike momentum axis once the additional scalar coupling is taken 
into account.^ Our results for the pole masses obtained in these models are given in Table HI 
For the model interaction Eq. ()28|) we agree with the estimate for the singularity closest 



® We liave even arbitrarily changed ce{0) from its value 2.97 in these fits. Dynamical chiral symmetry 
brealcing occurs for a(0) > Ofcrit with Ofcrit being slightly below one. For a(0) in the range Ofcrit < a(0) < 10 
we found no evidence for positivity violation when the CP vertex is used.^ 

^ Note that a similar result has been found in the model study of Ref. |43| where a Stingl-type gluon 
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FIG. 6: The dimensionless Schwinger functions msingAt,(t) (left) and As(t) (right) as function of 
msing t, where rngmg is the "pole mass" as determined by the exponential decay of the Schwinger 
function for different current quark masses m^, renormalized at /x = 10 GeV. 



to = given in Ref. [2^ based on a Taylor series expansion of the quark propagator 
functions, confirming that we can indeed extract the location of the first singularity via 
the Schwinger functions. Finally, we checked the truncation scheme of Ref. j75[ where a 
model interaction with an infrared finite coupling has been employed together with a bare 
quark-gluon vertex. In this case we also found a pair of complex conjugate poles as could 
be expected. 

Another interesting property of expression ()25|) is its insensitivity to explicit chiral sym- 
metry breaking, i. e. a current quark mass. The contributions from current quark masses to 
the function B{p'^) are almost momentum independent and therefore cancel quite accurately 
in Eq. (j25p . The Schwinger functions become steeper with increasing quark mass, but show 
no signs of positivity violation, even for current quark masses as large as a few GeV. For a 
detailed comparison of the mass dependence of the Schwinger functions ^sif) and A.y(t), we 
scale At,(t) by the pole mass, mging (extracted from the exponential decay of As^^(t)), and 
plot As(t) and rngmg A^(t) as function of the dimensionless variable mging^ in Fig. El This 
reveals that the only mass dependence is in the curvature of As(t) at small rngmg with 
increasing current quark mass the amount of curvature decreases. 

How can we understand this curvature that is present in Agit) but not in A^(t)? A 
possible origin could be the fact that the function (Js{p^) drops off like in the chiral 
limit while o"„(g^) decreases as As can be seen from Eq. (jH)), a single real pole on the 



propagator model has been employed in the quark DSE together with a quark-gluon vertex consisting 
only of the Abelian Ball-Chiu and Curtis-Pennington type structures |55|, |56j . In this study the absence 
of complex singularities in the quark propagator has been attributed to the vanishing of the employed 
model gluon propagator at zero momentum. This interpretation seemed to be supported by a study 
using the same propagator and a bare vertex which finds also real poles 0|. However, the present study 
clearly demonstrates that for a sufficiently strong interaction the crucial reason for this absence of complex 
singularities lies in the quark-gluon vertex. 
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negative momentum axis results in a pure, exponential decay of the corresponding Schwinger 
function. However, the Schwinger function of a propagator with two poles is 

oo 

- [ dpcositp) , ^ , , ^ = — — ^ ( -e"'"* - ^e"^*^ , (39) 



and for A somewhat larger than m, this could lead to the observed curvature at small t. This, 
in combination with the fact that this curvature tends to decrease with increasing current 
quark mass, suggests that this curvature is related to the l/p^ fall off (up to logarithmic 
corrections) of M(p2) in the chiral limit. However, there are other mechanisms that could 
generate such curvature as we will discuss in more detail in the next section. 

Comparing the two panels of Fig. El we also see that As(t) approaches msingA^(t) from 
below for all values of the current quark mass. In other words, we find that (within nu- 
merical accuracy) msingA^(t) > As{t) for all t. Based on the constraint for the spectral 
decomposition, Eq. (fT^ . this is what one would expect for a propagator describing a Dirac 
field with asymptotic states. Thus, within this approach there are no signals of positivity 
violation in the non-perturbative quark propagator. 

Considering these findings, we state the second major result of this work: the presence 
of a scalar quark-gluon coupling of sufficient strength leads to a positive definite quark prop- 
agator with a singularity on the timelike real momentum axis. As our quark-gluon vertex 
has been constructed as an ansatz, we do not have model independent information on the 
relative strength of the different tensor structures in the true quark-gluon vertex. Our as- 
sumption has been that all non-Abelian corrections can be accounted for by an overall factor 
multiplying an Abelian construction for the tensor structure of the vertex (see Eq. (H^). 
This factorization assumption has been tested in a recent investigation of the quark-gluon 
vertex in quenched lattice QCD and was found to be only valid at a qualitative level [t^. 
However, as yet no definite statements can be extracted from the lattice calculations as 
they are only performed in two special kinematical situations, whereas in our calculations 
the vertex is probed over the whole range of momenta. Further investigations are neces- 
sary to determine the relative strength of the various components of the vertex in a model 
independent manner. 

In QED4 however, we encounter a somewhat different situation. The vertex construction 
is more constrained than in QCD as the longitudinal part of the CP vertex, the Ball-Chiu 
vertex [H^, is exact and the relative strengths of the three longitudinal Dirac structures 
in the vertex are uniquely determined by the Ward identity, Eq. (j24p . The results for 
the fermion propagator in quenched approximation (a(g^) = a, constant) in the chirally 
broken phase of quenched QED are very similar to those of QCD. Again, we find a fermion 
propagator that satisfies positivity as long as it is calculated with a vertex obtained from 
the Ward identity but violates positivity if a bare vertex is used. The Schwinger functions 
are shown in Fig. [7| and the deduced (complex) pole masses are included in Table HI Of 
course, it remains possible that the transverse parts of the exact vertex conspire to lead to 
positivity violation again. However, this is unlikely, in particular in QED where one has no 
confinement. 
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FIG. 7: Results for two different vertex constructions in QED4. For ease of comparison we 
employed two different values for the coupling, i.e. a = 1.2 in the case of the bare vertex and 
a = 1.06 for the case of the CP vertex. 

IV. ANALYTIC PROPERTIES OF THE QUARK PROPAGATOR FROM PA- 
RAMETERIZATIONS 

In this section we explore the possible analytic structure of the quark propagator in more 
detail. Here we also consider the available lattice data for the quark propagator and investi- 
gate whether it is possible to obtain information on the analytic structure of the propagator 
by fitting this data, the DSE solutions, and the corresponding Schwinger functions with 
different parameterizations of pole locations and/or branch cuts. The singularity on the 
real momentum axis may be accompanied by additional real singularities at larger mass 
scaled or by complex conjugate singularities with a larger real part of the mass, or it may 
be the starting point of a branch cut on the negative real momentum axis. In the next two 
subsections we explore these possibilities. 

A. Meromorphic parameterizations 

The most rigorous constraint on the non-perturbative quark propagator is that it must 
reduce to a free fermion propagator at large momenta because of asymptotic freedom. This 

entails that the propagator functions, Cs^vip^) — ^ in all directions of the complex p^- 
plane (t^]. Additionally, the theory of complex functions tells us that if o"t,(p^) and <7s{p'^) 
are not constant, they cannot be analytic over the whole complex plane: non-constant, 
entire functions which are analytic at all finite points in the complex plane are already 
excluded by the asymptotic properties of the propagator functions. From the truncated set 
of DSEs explored in the previous section, we found the dominant (in terms of the Schwinger 
function) structure to be either a singularity on the negative real axis or a pair of complex 
conjugate singularities in the left half of the complex p^-plane. In both scenarios the poles 
are accompanied by additional undetermined structures which are responsible for the small 
time behavior of As(t). Guided by these results we first consider parameterizations of the 
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renormalized quark propagator using the meromorphic form 

np 



np / s 

^ v*/^ + «j + + «i - *i ) 



(40) 



with Tip pairs of complex conjugate poles located at ± ihj with residues r^. This form 
includes the possibility of complex conjugate as well as purely real poles, but enforces neither 
of these from the outset. Similar simple parameterizations have been considered in Refs. j3^ . 

In the following, we use physical constraints as well as lattice data to fix the position of 
the various singularities. The only practical restriction on this procedure is in the number 
of parameters that can be pinned down. As further simplifications, we assume that the 
residues, r^, of these poles are real (although this is not a strict requirement) and only 
consider the chiral limit. 

For the propagator functions, o"^(p^) and cr„(p^), the form Eq. ()40|) simplifies to 



np 



In terms of these quantities, we can construct the usual renormalization point independent 
mass function M(p^) = 0"s(p^)/cT„(p^) and the wave-function renormalization Zf{p'^) = (p^ + 
M^(p^))(Tt,(p^). In order to make contact with lattice data (where the finite lattice spacing 
leads to a maximum possible momentum), we renormalize at /x^ = 16 GeV^. 

There are various restrictions we can impose on the parameters rj, aj and bj in the mero- 
morphic form, Eq. (j4(J|) . These arise from its mathematical properties, from experimental 
observables and from recent lattice data. Asymptotic freedom requires that quarks behave 
like free particles at large momenta. Consideration of the large momentum limit of cr^(p^) 
implies that 



np 

E^^ = o- (43) 



Since we are working in the chiral limit, the mass function, M{p'^), must vanish for large 
spacelike real momenta. This entails that^° 



np 



J]r,a, =0. (44) 

Furthermore, M{p'^ +oo) must be real and approach zero from above. 
Asymptotically, the chiral limit mass function behaves as 

M(r)'^~^^-=^L_, (45) 



Nc 2 



^ -QCD 



^'^ If we move away from the chiral hmit, the right hand side of Eq. 14411 is replaced by the renormalized 
current mass. 
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where {qq) is the renormahzation-point-invariant chiral condensate. Although the loga- 
rithmic behavior of Eq. ()45p cannot be reproduced by these simple meromorphic fits, the 
logarithm is a slowly varying function and we estimate the condensate by fitting the mass 
function with Eq. over the range G (10^, 10^) using the Aqqd = 0.5 GeV and the 
appropriate 1-loop value of ■jm = 12/33 for Nf = 0. We then insist that this conden- 
sate extracted from our meromorphic propagator agrees with the phenomenological value 
{qq) = -[0.275(75) GeVf. 

In order to be phenomenologically apphcable, the propagator should reproduce the pion 
decay constant to a reasonable accuracy. To calculate this, we employ the approxima- 
tion [53], 

A2 



^ Z2-^ / dp^ 



2^ 



47r2 7 ^ ^ Zf{p^ 



p^ [ da^{p'^) 2n / 2\d(^s{p^ 







2 \ dp"^ dp"^ 



(46) 



which incorporates only the effects of the leading Dirac structure of the pion Bethe-Salpeter 
amplitude in the chiral limit. From a comparison of the relative sizes of the pion Bethe- 
Salpeter amplitudes in model calculations 0,1^, one concludes that this approximation 
should lead to an underestimation of /^r by 10-20 %}^ In our meromorphic fits we therefore 
demand that Eq. (gg) gives ~ 0.08(3) GeV. 

The Landau gauge quark propagator has been investigated on the lattice by a number 
of different groups using mean- field- and non-perturbatively- improved clover actions [t^, 
the Kogut-Susskind action j?^, the overlap formalism [t^ and the Asqtad quark action 
[t^I . The data sets obtained in the latter two formulations have the smallest error bars and 
are therefore employed in what follows. Their mass functions and wave-function renormal- 
izations have already been shown in Fig. |3] The mass function data from the lattice have 
been quadratically extrapolated to the chiral limit, whereas the mass dependence of 



Z^{p^) is very mild so no extrapolation has been performed. While the simple extrapolation 
procedure that has been employed may lead to sizable errors it will prove sufficient for 
our purposes. 

Unfortunately all of the lattice studies to date make use of the quenched approximation. 
Removing all internal quark loops is a potentially drastic modification of the theory. It 
destroys the unitarity of the S-matrix, however it is often assumed that these violations of 
unitarity are small. Strictly speaking, it is nonsensical to discuss the concept of positivity 
in such a situation and the lattice data discussed above cannot be relied on to provide any 
guidance in studying positivity of the quark propagator. However, from our experiences with 
the DSE studies of the previous section, one may expect that quenching will not qualitatively 
change the momentum dependence of the propagator (see Fig. 0} . Additionally, the lattice 
data apparently still contain large finite volume effects (especially in the wave-function 
renormalization) [s^, and do not precisely constrain the asymptotic (jP' +oo) behavior of 
the propagator . For these reasons we do not directly fit the lattice data (though a posteriori 
fits to it return very similar parameters to those we find below), but merely extract its 
three qualitative infrared features. Thus we assume that the zero momentum values of the 

7f 

^0) 



mass function and wave-function renormalization, Mq and Zq , and an approximate width of 



One also knows from chiral perturbation theory that the chiral limit pion decay constant is somewhat less 
than the physical value of 93 MeV. 
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n 


ai [GeV] 


hi [GeV] 


r2 


aa [GeV] 


62 [GeV] 


rs 03 [GeV] 


3R 
2CC 

IR+ICC 


0.365(15) 
0.360(22) 
0.354(15) 


0.341(25) 
0.351(69) 
0.377(64) 


0.08(5) 


1.2(8) 
0.140(*) 
0.146(*) 


-1.31(12) 
-0.899(*) 
-0.91(*) 


0.463(75) 
0.45(7) 


-1.06(*) -1.40(*) 



TABLE II: Best fit parameters of tlie tliree meromorphic forms: three real poles (3R), two pairs of 
complex conjugate poles (2CC) and one real pole and one pair of complex complex conjugate poles 
(IR+ICC). The parameters whose errors are replaced by an asterisk, are completely determined 
in terms of the other parameters through Eqs. ()43() and 1)441) . In order to reproduce the results 
presented here, one should use the values that follow from Eqs. ()43j) and ()44|) for those constrained 
parameters. 

the region of large dynamical mass generation, (defined by M{ijj\) = Mo/2), are robust 
against the effects of quenching (within substantial errors). With this in mind, we require 
that our parametric fits are in reasonable agreement with the extracted values of Mq, Zq 
and Ul. That is: 

Mo = 0.35(10) GeV , Z/^ = 0.6(2) , ujl = 0.7(2) GeV . (47) 

Note that {qq), fn and these three parameters are obviously not entirely unrelated. 

Given the number of independent constraints we can impose, we can reasonably expect 
to be able to determine only five or six parameters. This implies np < 3 in Eq. (00)). We 
find that three paradigmatic cases satisfy the requirements of Eqs. (jlSjl -ipT j) : three purely 
real poles (denoted, 3R), two pairs of complex poles (2CC), and a real pole plus a pair of 
complex conjugate poles (IR+ICC). In order to construct the best fits for each of these 
forms, we first impose the simple constraints of Eqs. (pSj) and ()44|) to reduce the number of 
parameters to be varied. Then for each parameterization we randomly sample the available 
parameter space, constructing a large ensemble of parameter sets that satisfy the full set of 
constraints. The best fit parameters and their errors are finally calculated as the mean and 
standard deviation of the parameters in this ensemble. 

The simplest possible parameterizations of a single real pole or a single pair of complex 
conjugate poles (rip = 1 in Eq. (HO))) cannot satisfy the required constraints. Specifically, 
enforcing the perturbative asymptotic behavior (Eqs. PHjl and ()44|) ) makes it impossible to 
satisfy any of the other requirements described above. Similarly, for two real poles (rip = 2, 
61 = 62 = 0), the restrictions on the infrared properties (/^r, Mq and ujl) are incompatible 
with a realistic quark condensate. 

As mentioned above, a satisfactory realization of the requirements of Eqs. ()43j) - ()47p is pos- 
sible in the case of three real poles (rip = 3 and 61 = 62 = ^3 = 0). The best fit parameters we 
obtain are shown in Table ITTl and related quantities that they result in are given in Table ITTTl 
Although the propagator functions crs,^(p^) have poles at ~ —0.2 GeV^, they exactly 
cancel in the combinations M(p^) and Z-^(p^). However the functions M(j9^) and Z-^(p^) do 
have poles further in the timelike region, the first one occurring at ~ —0.75 GeV^. Also 
the zeros of Z^{p'^) on the real axis may be problematic as they will necessarily produce 
singularities in the CP construction of the quark-gluon vertex, c.f. Eq. ((221) • 

In the case of two pairs of complex conjugate poles (rip = 2), the best fit parameters and 
calculated quantities are again given in Tables ITTl and UTTl Both M(p^) and Z-^(p^) exhibit 
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Mo [GeV] 


7^ 


[GeV] 


U [GeV] 


[GeV] 


3R 


0.29(1) 


0.55(7) 


0.79(4) 


0.071(3) 


0.3(2) 


2CC 


0.33(11) 


0.57(12) 


0.69(27) 


0.070(31) 


0.3(3) 


IR+ICC 


0.31(7) 


0.52(7) 


0.72(25) 


0.068(23) 


0.3(2) 



TABLE III: Values for the various constrained quantities for the three parameterizations of Ta- 
ble ^ Errors are solely due to uncertainties in the parameterizations and do not include any 
additional systematic errors. 
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FIG. 8: Propagator functions for the fit using one real pole and one pair of complex conjugate 
poles. 

unexpected behavior around ~ —0.12 GeV^, where they have a very sharp pole and a zero 
on the real axis. This arises because (Jsip^) and <J^{p^) have zeros at very slightly differing 
momenta {p^ = -0.127 GeV^ vs -0.117 GeV^) and it may be somewhat troublesome. This 
behavior, as well as the small imaginary part of the location of the first pair of poles, suggests 
forcing the first pair of poles to collapse to one real pole (np = 2, 6i = 0). 

Redoing the fits with one real pole and one pair of complex conjugate poles, we come 
up with very similar parameters to the 2CC parameterization, as listed in Table |n] The 
corresponding propagator functions are shown in Fig. |H1 With this parameterization, the 
strange behavior of M(p^) and Z^{p'^) disappears and Z^{p'^) only has complex conjugate 
poles and zeros (Z-^(p2 = -0.41 ± 0.48z GeV^) = 0, Z^lp"^ = -0.55 ± 0.69i GeV^) ^ 
oo) so the longitudinal part of the quark-gluon vertex, Eq. (j23|) will not have particle-like 
singularities '81|. This parameterization also contains one parameter less than the others. 
Therefore we consider this to be the preferred form of the meromorphic parameterizations 
investigated here. 

In comparing the three sets of parameterizations, it is worth remarking that the location 
of the (real part of the) first pole and its residue are extremely robust. The obtained value 
for this constituent quark mass, m = 377(64) MeV for our best fit, is also in good agreement 
with a value extracted from lattice simulations of the quark propagator using a tree-level 
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■ ■ Fit: 2 CC poles 




p[GeV] 



FIG. 9: The best fit curves for the meromorphic parameterizations compared to the lattice data. 
For the details of the parameterization of crg^y^p^) by a form with a branch cut, fitted to the Asqtad 
data, see the next subsection. 

Symanzik improved action, m = 342(13) MeV js^l- However, the constraints on the other 
features in the fits are less precise, especially in the case of three real poles. In Fig. El we 
compare the parameterizations given in Table ITTl to the lattice data; overall, the agreement is 
quite acceptable. Note that the meromorphic fits have relatively low values of Z^; this may 
change once finite volume effects are reduced in the lattice data. Also, each parameterization 
has a somewhat low value of f^^ in the chiral limit. This can be attributed on the one hand 
to the approximation leading to Eq. and on the other hand to the approximations on 
the lattice: the chiral extrapolation as well as the omission of dynamical quarks might lead 
to an underestimation of in the lattice data jl3l |. 

Having determined the best parameters for three different forms of our fit functions, we 
now examine the Fourier transforms of the momentum space propagator functions as^vip^)- 
Specifically, we attempt to determine whether the sub-dominant behavior of the various 
parameterizations can be determined from the Schwinger function, and, if so, apply this to 
the DSE solutions of Sec. 11111 

Using the identity 

oo 

d^^rr^ = ^ ^M-cy) [y e R, arg(c2) ^ vr] , (48) 



we can directly calculate the Schwinger functions from our parameterizations (Eqs. ()4H1 and 

my- 

np 

A,{t) = ^s^n(ai)rie-l"»l*cos(6it), (49) 

i=l 

^v{t) = 2J 2 I h2 cos{bit) - bi sm{bit)) . (50) 

i=l + "i 



27 




t [GeV"'] t [GeV"'] 



FIG. 10: The Fourier-transform of both o"s(p^) (left) and C7^(p^) (right) for the optimal parame- 
terizations with multiple poles. In the diagram on the left hand side we have included our DSE 
results. 

For all parameterizations, the term with the smallest mass parameter a, will dominate for 
large t. In Fig. ^1 we display the analytic Fourier transforms of the parameterized scalar 
and vector propagator functions, Eqs. and (jnOj)- For comparison, we also include our 
DSE result for As(t) employing the CP vertex. Note the qualitative difference between 
the parameterization with two complex conjugate poles and the other two. Whereas the 
latter show no sign of positivity violation, in the 2CC parameterization we clearly see zero 
crossings of the Schwinger functions, both in As{t) and in Ay{t) (even a small imaginary 
component in the complex conjugate masses is detectable provided the Fourier transform 
can be calculated accurately to large enough t). Note that As{t) calculated from the mero- 
morphic parameterizations shows a similar amount of small t curvature to the DSE result, 
but Ay{t) is linear in this region. Thus multiple poles as explored here could explain the 
small t behavior observed in the DSE Schwinger functions. 

We also use these analytic Fourier transforms to test our numerical Fourier transform, 
finding that it reproduces the analytic results down to As,„(t) ~ 10~^ where we begin to 
run into accuracy problems. However, the numerical routine we employ is clearly able to 
distinguish between a dominant real pole and dominant complex conjugate poles. This gives 
us further confidence that our results from the DSE solutions in the previous section are not 
numerical artefacts. 

B. Parameterizations with branch cuts 

As mentioned above, there is evidence that the Schwinger function As{t) is convex (with 
sizable curvature) at small t. On the other hand, the Schwinger function Ajj{t) as obtained 
from the DSE solution shows no such curvature. This difference could be accommodated 
within the simple meromorphic fits of the previous subsection. However, this is certainly 
not the only possible mechanism leading to such a difference, and here we explore the 
consequences of allowing for singularities with branch cuts. As can be seen from Fig. ^ the 
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curvature of As(t) depends on the current quark mass, so we also consider the effects of 
exphcit chiral symmetry breaking. 

Our motivation for investigating such parameterizations arises from considering the DSE 
for the quark propagator, Eq. (|27|). If the combination a{k'^)/k'^ is non-analytic at fc^ = 
(g — pY = (in other words, if a{0) 7^ 0), the integration path necessarily passes through 
the external point p. Thus, in order to evaluate the quark propagator at arbitrary complex 
momenta, one has to deform the integration contour in the DSE and solve the DSE along this 
deformed integration path. As long as there are no singularities in the other factors of the 
integrand (i.e. in S{q) and V^'^^^q, k)), this can in principle be done unambiguously (though 
it is numerically a nontrivial task). However, if we want to evaluate the integral for a value 
p at which the propagator, S{p), has a singularity, we are forced by the analytic structure of 
a{k'^)/k'^ to include this value of p in the integration contour for d'^q. Thus, we have a pinch 
singularity at this point coming from S{q) and a{k'^)/k^] this generally leads to a branch- 
cut, as is shown in more detail in Appendix B. We also note that the asymptotic form of 
the quark propagator has perturbatively calculable logarithmic contributions. Considering 
these points, we would expect that the singularities in cTg^vip'^) are branch points rather than 
simple poles. Thus we next attempt to parameterize the quark propagator by functions with 
branch cuts using the parameterization of the strong running coupling, Eq. (jH^ . that has 
proven helpful in understanding the analytic structure of the gluon propagator. 

As a first try, we shall fit the inverse propagator functions M(p^) and {p^) as obtained 
from the quark DSE with the CP vertex. Given the close agreement of the DSE solutions and 
the lattice quark propagator seen in Fig. |21 fitting the DSE solution will result in similar 
physical constraints to those of the previous subsection. The leading-order perturbative 
behavior is known, and we allow for one additional sub-leading term, that is to be fitted to 
the DSE solution. Furthermore, we want the parameterization to have a branch cut along 
the negative real axis starting at p^ = —m^i^g- Thus we are lead to fit the DSE solutions 
with 

ZHp ) = ^1 — + -r^^T^.) . (51) 

M{p') = C,^s, + , 2^ 2 \ A2,2 + Cc,m a{p' + ml^^y- . (52) 

P + "^sing + ^ [P + "^sing + ^ ) 

The parameters Cdcsb and Ccqm are related to the chiral condensate and the renormalized 
current quark mass, respectively: 

- {m) = c,c., [^) ^ , (53) 

r I 

HpV^hcD) 



m 



The renormalization constant Z2 is determined by the renormalization condition {^x^) = 1, 
msing follows from the exponential decay of the Schwinger functions, and we take A to be 
equal to Aqcd in the running coupling, a{x), for which we use Eq. (IH^ . The remaining free 
parameters in this fit, C2 and C4, are fitted to the numerical solution of the DSE and Aqcd 
is also varied to improve this fit. 

The results are shown in Fig. ^2 for several different current quark masses representative 
of masses up to that of the bottom quark. The fitted parameters are given in the first section 



29 




FIG. 11: The functions (p^) (left) and M(p^) (right) obtained from parameterizations fitted 
to the numerical DSE solutions for different current quark masses. The solid curves are the fits 
(Eqs. (|5T|) and (|52|) ) and the dotted, dashed, and various dot-dashed curve are the DSE solutions 
as shown in the legend (the curves for Z-^(p^) in the chiral limit and for = 5 MeV are almost 
indistinguishable). The fit parameters are given in Table ITVl 
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TABLE IV: Parameters for the various branch cut fits (Eqs. ([5T|) - (|60|) ) to the Nf = quark DSE 
solutions for different masses. A parameter set obtained by fitting Eqs. ([55 |l -(|57 j) to the Asqtad 
lattice data is also included. 
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FIG. 12: The Schwinger functions A<j(t) (left) and Ay{t) (right) for the three parameterizations 
with branch cuts, in the chiral Umit. For comparison, we have also included our DSE results. 



of Table HVl With only a few parameters, the fits represent the DSE solutions very well over 
the entire Euclidean region. The fitted values of Ccqm are all reasonably close to the current 
quark masses that were used as input in the DSEs (small deviations are due to sub-leading 
effects) and the (fitted) chiral condensate is acceptable: —{qq) = (290 MeV)^. 

Despite the fact that these parameterizations fit Z^ijP) and M{p^) so well, the cor- 
responding Schwinger functions do not fit the Schwinger functions of the DSE solutions. 
Clearly, the zeros of + M'^ijP) (which determine the poles of as^viv^)) will in general not 
occur on the negative axis when Eq. (j52|) is used for M{p'^). Indeed, the dominant singu- 
larities of the propagator functions as^vij)^) calculated from the parameterizations of (j)^) 
and M{p'^) are a pair of complex conjugate singularities, and the corresponding Schwinger 
functions clearly show oscillations^^ (see Fig.ll2|). Extensive "fine-tuning" of the fitting form 
and/or the parameters is required in order for p^ + M'^{p^) to have its first zero at the pole 
mass deduced from the Schwinger function of the DSE solution. 

As a second alternative, we can directly parameterize as,v{p'^), and fit these to the nu- 
merical DSE solutions. Again, we want to reproduce the leading logarithmic corrections to 
Z-^(p^) and M(p^), which can be achieved by using the forms 



-Bchiral(P^) 



^sip') 
<^v{p'^) 



c, 



dcsb 



+ 



A2)2 ' 



5chiral(p^) , Ccqm a{p 



m. 



P + TTlgjjjg 
1 

p2 + TTigj^g 



+ 



sing; 



p^ + m. 



(55) 
(56) 



smg 



1 - 



27r 



C i?chiral(p') ) . (57) 



This form has mass-like singularities in (Ts^vip"^) at p"^ 



— "^oino^ from which branch cuts 



^"^ For the heavier quarks, these oscillations are numerically difficult to detect because the Fourier transform 
falls off very rapidly with t. 
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p^[GeV^] p^[GeV^] 

FIG. 13: The functions (jp') (left) and M{p'^) (right) obtained from parameterizations of 
Us^vip^)-, Eqs. (|53)) -(|57 |) . fitted to the numerical DSE solutions for different current quark masses. 
The solid curves are the fits, with the fit parameters given in Table HVl 



extend to = — oo. Away from the real axis, Cs^vip"^) have no singularities (though there 
is a second singularity at p"^ = — wigj^g — A^). Furthermore, this parameterization ensures 
the correct asymptotic behavior, both for crs,v{p'^) and for the quark functions M{p'^) and 

ijp'). The main disadvantage of fitting o's^yijp') is that the analytic structure of ijp') and 
M{p^), and of A{p^) and B{p^) will now become non-trivial. Again, a delicate fine-tuning 
is required to obtain a good fit for both a,,„(p2) and for Zf{p^), M{j?), A{p^) and B{j)^). 

The parameters play a similar role to those in the previous parameterization, with the 
exception of C which is determined by requiring that Z^ijP') is finite at the mass pole. 
The other parameters are fixed by fitting Z^i^p"^), M^p"^), and the Schwinger functions. For 
moderately small current quark masses, we can obtain reasonably good fits, as can be seen 
in Fig. Uni with the corresponding parameters listed in Table HV1 We can also fit the Asqtad 
lattice data quite well with this parameterization, as shown in Fig. El For current quark 
masses larger than a few hundred MeV, the wave function renormalization can no longer be 
fitted with this relatively simple form. This is most likely related to the substantial increase 
in the constant C2 for heavy quarks when fitting Z^ijP') directly (see Table HV]) . 

The functions v4(p^), B{p'^), MijP') and Z^{p'^) have a singularity at p^ = — wigj^g where a 
branch cut along the negative real axis starts, and another singularity further in the timelike 
region at = — rngj^g — A^. In addition, Mij?) and Z^ {p^) have a pair of complex conjugate 
poles located at the zeros of a^ijP'), and v4(p^) and -B(p^) have two pairs of complex conjugate 
poles at zeros of p^ crtip^) + o"s (P^)- 

The Schwinger functions, /S.s,vit)-i are reproduced very well, see Fig. [T21 Notice that the 
parameterizations of Zf{p^) and M{p'^) fit the DSE solutions for Z^{p'^) and M{p'^) better 
than these parameterizations of ag^vip^)-, whilst the latter parameterizations are obviously 
better fits of the Schwinger functions corresponding to those same DSE solutions. Thus we 
are warned that even an almost perfect fit for Euclidean momenta does not guarantee a 
good fit of its Fourier transform, let alone a good representation of the function in the entire 
complex plane. 
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FIG. 14: The functions {jp') (left) and M{p^) (right) obtained from parameterizations of A^p'^) 
and B(p'^), Eqs. (|58() -()6U |) . fitted to the numerical DSE solutions for different current quark masses. 
The solid curves are the fits, with the fit parameters given in Table HVl 



Finally, we construct a parameterization of the inverse quark propagator functions A{p'^) 
and Bij?), such that the propagator functions (Ts,v{p^) have pole-like singularities on the 
timelike p^ axis. For this purpose we use the parameterization 



with 

. (C^, ^ + g + C.,.. a(0)-) / (1 + + §) , (60) 

and Z2 determined by the renormalization condition A{fi'^) = 1. The results of these fits 
are shown in Fig. El with the corresponding parameters given in Table IIVI Though not as 
good as the direct fits of (jp') and MijP') (Eqs. (jSip and (j3^ ) that were made without 
taking into consideration the analytic structure of ag^vip^), these fits reproduce the DSE 
results within about 10 to 20% over a wide range of masses. By construction, the dressing 
functions again reduce to the perturbative forms in the ultraviolet region and the analytic 
structure is in agreement with the Schwinger functions corresponding to our DSE solutions. 
In the chiral limit, As(t) corresponding to these fits shows significant curvature at small t, 
as can be seen from Fig. for larger quark masses this curvature decreases. In contrast, 
A.t,(t) does not show this curvature in agreement with our DSE results. However, the actual 
analytic structure of (Js,v{p'^) is rather complicated. In addition to the singularity on the 
negative real axis at = —"nosing "where a branch cut starts, it also has a pair of complex 
conjugate poles at zeros of p"^ A^(p^) + B'^ijP'). 
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C. Generic features of the quark propagator 

The results of this section point strongly towards an analytical structure of the quark 
propagator with a dominant singularity on the real timelike axis. At present, the nature of 
this singularity cannot be determined with confidence. It could be a simple pole, in which 
case additional poles or other types of singularities further away from p^ = are needed 
in order to explain the observed behavior of the Schwinger function. However, given the 
structure of the quark DSE, it is more likely that this singularity is a branch point, and 
that there is a branch cut along the negative real axis starting there. Having only one 
branch point singularity on the negative real axis is (in principle) sufficient to reproduce 
the observed Schwinger functions. We have not yet been able to distinguish between these 
alternatives by numerical calculations of the Schwinger functions of the DSE solutions. 

Given the strong sensitivity of the Schwinger functions to the details of the propagator 
functions, and the fact that the dominant singularity is well into the timelike region, it is un- 
likely that the sub-dominant analytic structure of the quark propagator will be determined 
by Euclidean lattice simulations. The situation for the gluon propagator is quite different: 
there, the analytic structure is highly constrained by the behavior for — > 0"*". By approach- 
ing the singularity at = from the spacelike region, we can gain information about the 
nature of this singularity. In contrast, the first non-analytic point of the quark propagator 
is (most likely) at = —m^-^^g < 0. Thus, the behavior of M{p'^) and Zf{p'^)^ or as^vip^)-, for 
p^ — > 0+ does not reveal much about the analytic properties of the propagator. We cannot 
approach the singularity without accessing the timelike region explicitly. In addition to this 
dominant singularity on (or very close to) the negative p^ axis, the propagator may have 
other sub-leading singularities further away from p^ = 0. Within the DSE framework one 
would have to solve the quark DSE over a suitable region of the complex momentum plane 
to decide questions about the nature of the dominant singularity and about the existence of 
sub-leading singularities further from = 0. However, this is numerically very demanding 
and not within the scope of the present investigation. As we have seen in the previous 
section, the results could also be strongly infiuenced by the truncation of the DSEs. 

In all of the parameterizations of the preceding subsections a robust feature appears: the 
leading singularity is on (or very near to) the real axis. The scale at which this singularity 
occurs is somewhat dependent on the constraints used in the fits; the lattice data suggest a 
scale of 350 to 390 MeV, whilst the DSE solutions prefer a slightly larger value ~ 500 MeV. 
Despite this slight variation, these numerical values are somewhat intriguing and hint at a 
possible interpretation in terms of a constituent quark mass. 

Fits like the examples presented here might be usefully applied in hadron phenomenology 
pending a more conclusive investigation of the analytic properties of the quark propagators. 
However, one should treat these parameterizations with care and keep in mind that neither 
the fitting forms nor the parameters are unique. The same Euclidean data, from lattice or 
DSE calculations, can be fitted quite well with different parameterizations having distinct 
analytic properties. The only robust feature concerning the analytic structure is that the 
dominant singularity, as probed by the Schwinger function, is on (or very near) the real 
timelike axis. 
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V. SUMMARY AND CONCLUSION 



In this work we have investigated the analytic structure of the propagators of Landau 
gauge QCD in the timehke momentum region using Schwinger functions, and employing 
various analytic parameterizations. We summarize the main results below. 

Both lattice simulations and Dyson-Schwinger equation calculations suggest that the 
gluon propagator is finite or even vanishes in the infrared. The latter behavior necessarily 
leads to violations of reflection positivity, a sufficient (but not necessary) condition for gluon 
confinement. Our numerical analysis of the Schwinger functions calculated from the DSE 
solutions confirms this behavior, finding clear evidence of such positivity violations in the 



gluon propagator in accordance with previous results p, [111, |21[. The gluon Schwinger 
function possesses one zero at t ~ 1 fm, marking the length scale above which sizable 
negative norm contributions appear. We explore the detailed analytic structure of the gluon 
propagator in the timelike {p^ < 0) momentum region by constructing parameterizations 
that fit both the spacelike momentum behavior of the lattice calculations and DSE solutions, 
and the corresponding Schwinger function. These parameterizations incorporate the power- 
law infrared behavior determined analytically from the coupled ghost and gluon DSEs, and 
the perturbatively calculable ultraviolet logarithmic behavior. The crucial feature of these 
parameterizations is the presence of a branch cut on the timelike momentum half-axis which 
produces the observed positivity violations. These simple parameterizations depend on 
(effectively) only one parameter, the scale Aqcd- 

In exploring the analytic structure of the quark propagator using the same Schwinger 
function methods, we have found an unexpected sensitivity of this structure to the trunca- 
tion of the quark DSE. Gauge symmetry (or, more precisely, the relevant Slavnov-Taylor 
identity) requires the presence of a scalar coupling in the non-perturbative quark-gluon ver- 
tex. This coupling is only present if chiral symmetry is broken dynamically, in which case 
it leads to a self-consistent enhancement of the effective quark-gluon interaction. This term 
is important in at least two ways. First, it is required in order that the solution of the 
coupled quark, gluon, and ghost DSEs generates enough dynamical chiral symmetry break- 
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ing in the quark propagator to give a phenomenologically acceptable quark condensate 
More importantly for our investigations, the scalar coupling leads to qualitative changes 
in the analytic structure of the quark propagator, and hence has significant consequences 
for the questions of (non-)positivity and the manifestation of confinement. When this term 
is omitted (as in the commonly used rainbow truncation), positivity violations consistent 
with complex conjugate singularities in the quark propagator (as found in previous studies 
[13L 31, 113, 113) are unambiguously observed. However, when the gauge-mandated scalar 



coupling is included, no such evidence of positivity violation is found and the dominant 
analytic structure appears to be a singularity on the real, timelike (p^ < 0) axis. Whilst the 
absence of positivity violations says nothing about quark confinement (positivity violation 
is a sufficient but not necessary condition), it does tell us that confinement is probably not 
manifest at the level of the propagator. We also see similar behavior in (quenched) QED4. 
Here, a positive definite propagator is desirable as the electron is an observable particle. 

Finally in Sec. IIVI we have attempted to probe deeper into the analytic structure of the 
quark propagator. We have constructed various parameterizations and used lattice data, 
DSE solutions, and other general properties to constrain them. Whilst an infinite variety of 
functional forms (we have investigated only a few that come easily to mind - constructed 
from real or complex conjugate poles, and branch cuts on the timelike momentum axis) 
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would be capable of satisfying our constraint criteria, one robust feature emerges from our 
analysis: the dominant {i.e., closest to = 0) analytic structure occurs on (or very near 
to) the real, timelike half-axis. The scale of this mass-like singularity, as suggested by 
meromorphic parameterizations of the lattice data, is 350 to 390 MeV. The DSE solution 
indicates a scale of about 500 MeV. An accurate determination of the precise nature of this 
singularity and additional sub-dominant contributions awaits future improvements. 
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APPENDIX A: FURTHER DETAILS OF THE DSE TRUNCATION 

Since the quark-gluon vertex appears to play a significant role in the analytic properties 
of the quark propagator, we have explicitly explored the effects of this part of our truncation 
scheme in the main body of the text. Here we examine the various truncations we use in the 
Yang-Mills sector and investigate reasonable modifications to highlight the extent of both 
the truncation dependence and the truncation independence of our results. 

In this work, we have used the perturbative ghost-gluon vertex 

rf-*(g,p)=^g^. (Al) 

Lerche and von Smekal J§] have investigated a large class of possible structures for this 
Green's function. In particular their results show that such variation leads to infrared 
behaviour of the gluon and ghost propagators as in Eqs. (fTH|) and (fT^ with the exponent k 
in the range 0.5 < k, < 0.7. In this range the essential analytic properties that we find for 
the gluon propagator remain unaltered. 
For the three-gluon vertex, we use js^l 



P . , _ 1 G(g^)(^-/^-^°) Giiq - pf)i^~^/^-^^) ^ 

^p.a[q,P) ^^(^^^) Z(g2){i+a) Z((g-p)2)(i+^') ^-^^'^^ ^"^^^ 

where r^pJcr{q,p) is the perturbative form, 6 = —9Nc/{4:ANc — 8Nf) is the one-loop anomalous 
dimension of the ghost propagator, Zi is the three-gluon vertex renormalisation, and we fix 
a = b = 36. As discussed in the main text, this form is chosen to ensure the running coupling, 
Eq. fll6p . is renormalisation point independent and that the ghost and gluon propagators 
have the correct one-loop anomalous dimensions (these constraints are satisfied for arbitrary 
values of the parameters a and b). 
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Truncation 




It i 




It 11 




^ 1 


Knrn [MeVl 




Knrn fMeVl 

O LJ L J 


Full QCD, a = 36 




2.4(2.0) 


500(470) 


2.5 


510 


Quenched QCD, a 


= 36 


2.9 


550 


3.2 


550 


Quenched QCD, a 


= 26 


2.7 


410 


3.0 


400 


Quenched QCD, a 


= 46 


3.2 


560 


3.6 


550 



TABLE V: Gluon fit parameters for various truncations of the system of DSEs. The parenthesized 
results for the full QCD fit I are when the fit is optimized to reproduce fit the DSE Schwinger 
function; all other parameter sets are fitted to Z{p). 



Whilst no systematic study of this vertex truncation exists, in Ref. one of us has 
investigated some variation of the parameters a and h. There, the gluon propagator was 
found to vanish at — *■ independent of the choice of a and 6, and our conclusions 
about the analytic structure of the gluon propagator remain qualitatively unchanged. To 
investigate the dependence of the gluon fit parameters {ujiji and Is^qcd) on the three-gluon 
vertex truncation, we have made fits to the quenched DSE solutions using Eq. ljH^ for a 
number of choices of a and h. The resulting parameter sets are shown in Table |3 and are 
seen to vary by 20%. With certain choices of these parameters, one is able to closely match 
the quenched lattice data. However the original [a = h = 35) truncation gives results in 
adequate agreement (see Fig. 12)) with this data (which itself is not without error bars) and 
we use it exclusively in the main text. 

Finally, the truncation we employ neglects the effects of the four-gluon vertex. These 
effects are unknown up to now. Since the ghost loop is dominant in the infrared 4^ and 



the one-loop diagrams dominate in the perturbative, ultraviolet region, such effects are 
expected to be most important in the intermediate momentum regime {p^ ~ IGeV^). A 
two-parameter model for the corresponding two-loop diagrams in the coupled gluon-ghost 
DSEs has been explored in Ref. js^l ; under such variation the — > behaviour of the gluon 
and ghost propagator remains qualitatively the same as in our results. 



APPENDIX B: POLE OR BRANCH POINT IN THE QUARK PROPAGATOR? 

If the quark propagator has a non-analytic point at = — ?^sing where the propagator 
diverges, what kind of singularity can we expect? In order to answer this question, consider 
the generic integral that appears in the RHS of the quark DSE 

A 

Hp') = I d'q'^^^^^Ke{q\p\q-p) aU<f) ■ (Bl) 

Assume that the running coupling, a{k'^), does not go to zero like k"^ for k"^ —>■ 0. Further- 
more, assume that the kernel Kg{q'^,p'^, q-p) has no singularities. Thus the only singularities 
in the integrand are located at {q—pY = and at = — ^rigj^g (coming from the propagator 
function crs^(g^)). For Euclidean values of p^, we calculate this integral by performing the 
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x+iy 




FIG. 15: Location of the branch cuts (dotted curves) in K{q'^,p^) in the complex g'^-plane, for 
= X being real and for a complex value of = x + iy. Also shown are possible integration 
paths (solid curves) from to that do not cross the corresponding branch cuts. 



angular integral first, followed by the radial integral 

/(/) = j qW^sAQ')K{q',p'), (B2) 



with 

7r 

K{q ,p) = sm e de — Kg{q ,p ,qpcose) . (B3) 

J — 2qp cos u + p'' 



This angular integral is well-behaved for any Euclidean value of q^ and as long as the 
singularity in a{k'^)/k'^ is an integrable singularity, so let us assume from here on that this is 
the case. If we investigate the analytic properties of K{q'^,p'^) for arbitrary complex values 
of g^, while keeping real and positive, we discover that it has a branch cut that can be 
characterized by q^ = p"^ exp{i(j)) , with < (p < 2n. Notice that if we perform the radial 
integral along the positive real q^ axis, we do not cross this branch cut: the end-points in 
are not included. This is schematically depicted in Fig. IT^ 

Now consider the analytic continuation to complex values of p'^. Clearly, the branch cut 
in K{q'^,p'^) will move, as indicated in Fig. This means that the integration path in 
Eq. ()B2|) has to be deformed so as not to cross the (shifted) branch cut stemming from the 
angular integration, while keeping the end-points fixed. A possible integration path has been 
shown in Fig. though the actual integration path is of course not unique. The general 
rule for this deformed integration path is that it has to go through the point = p"^, since 
that is where there is an opening in the circular branch cut of K{q'^,p'^). This procedure 
leads to a well-defined and unambiguous analytic continuation of the Euclidean DSE, and 
can be implemented numerically p9L 

Following this procedure, one can now show that any singularity in o"s,t,(g^) leads (in 
general) to a branch point singularity in /(p^). This is shown in detail in Fig. ^Jfor the 
case in which o"s,„(g^) have singularities on the real time-like axis at = —m'^iT^g- We have 
drawn two distinct radial integration paths in order to calculate /(p^) for p^ = x — iy. One 
is obtained by continuously deforming the original integration path through the upper half 
of the complex g^-plane, crossing the negative real axis beyond g^ = —rnl^^^ (dashed curve). 
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FIG. 16: Two integration contours in the complex A;^-plane for the radial integral in /(p^) for 
= X — iy (left) and the resulting analytic structure for /(p^) in the complex p^-plane (right). 

and the other by deforming the integration path via the lower half of the complex plane (solid 
curve). Because of the combination of (i) the singularity in (Ts,«(g^) at = —m'^i^^ and (ii) 
the circular branch cut in K{q'^,p^), these two integration paths cannot be deformed into 
each other while keeping fixed. Therefore, the obtained values of /(p^) will (in general) 
be different, and I{p^) becomes a multi-valued function with a branch-point singularity at 

= — mgjng. The "natural" choice for the branch cut is along the negative real axis, as 
indicated by the wavy line in the right panel of Fig. [TBI 

Returning to the specific case of the quark propagator, we note that the RHS of the quark 
DSE contains an integral like /(p^), whereas the LHS is one of the inverse quark propagator 
functions v4(p^) or BijP'). If a^^si.'f) has a singular point k"^ = — "^smg' then /(p^) has a 
branch-point singularity at = —fn'^ij^g, and therefore A(p^) and B{p'^) will have a branch 
point at p^ = — "^gj^g. Thus, unless there are intricate cancellations, the singularity in 
o'v,s{p'^) of a self-consistent solution is a branch-point singularity and not a simple pole. 
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